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We study the question of prompt vs. delayed collapse in the head-on collision of two neutron stars. 
We show that the prompt formation of a black hole is possible, contrary to a conjecture of Shapiro 
which claims that collapse is delayed until after neutrino cooling. We discuss the insight provided 
by Shapiro's conjecture and its limitation. An understanding of the limitation of the conjecture is 
provided in terms of the many time scales involved in the problem. General relativistic simulations 
in the Einstein theory with the full set of Einstein equations coupled to the general relativistic 
hydrodynamic equations are carried out in our study. 
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Introduction. The study of the coalescence of neu- 
tron stars (NSs) is important for gravitational wave as- 
tronomy and high energy astronomy. However, at present 
we lack even a qualitative understanding of the process. 
One issue is the prompt vs. delayed collapse problem. 
While we expect that two 1.4 M NSs when merged 
will eventually collapse to form a black hole, the collapse 
could be delayed by fragmentation/mass shredding, an- 
gular momentum hang-up, and/or shock heating. The 
time scale of the collapse has important implications on 
the gravitational wave signals to be detected by LIGO 
jij. We focus on the issue of prompt vs. delayed collapse 
in this paper. 

The difficulty of getting an answer to this aspect of NS 
coalescence is very much the same as the full coalescence 
problem. Namely, we need to solve the full Einstein equa- 
tions coupled to the general relativistic hydrodynamic 
(GR-hydro) equations. 

Recently, Shapiro H put up an argument suggesting 
that one may be able to answer this question without nu- 
merical simulations, at least for the case of head-on col- 
lisions. The "Shapiro conjecture" goes as follows: Given 
the conditions: (I) that the two NSs are colliding head- 
on after falling in from infinity, and (II) the NSs are de- 
scribed by a polytropic equation of state (EOS) P = Kp T 
(with K a function of the entropy and the polytropic in- 
dex r remaining constant throughout the collision pro- 
cess), it is conjectured that no prompt collapse can occur 
for an arbitrary T and an arbitrary initial K. The ba- 
sic argument is that the potential energy when converted 
to thermal energy by shock heating is always enough to 
support the merged object, until neutrino cooling sets in. 

The argument based on conservation is appealing, and 
provides useful understanding for a range of the NS coa- 
lescence problems. However there is a major assumption 
for the argument to go through, namely, the collision 



process can be approximated by a quasi-equilibrium pro- 
cess, in two senses: (A) The coalescing matter can be 
described by one single EOS everywhere (K is a func- 
tion of time but not space), and (B) whether it collapses 
or not is determined by hydro-static equilibrium condi- 
tions, i.e., whether a stable equilibrium configuration ex- 
ists or not. This quasi-equilibrium assumption is not 
self-evident for the head-on collision of heavier NSs. It 
could happen that the coalesced object collapses before 
it can thermalize, or the collision process is so dynamic 
that even though a stable equilibrium state exists, it is 
not attained in the collapse process. The final outcome 
depends on the various time scales in the problem. 

Time Scale Considerations. We examine this as- 
sumption of "quasi-equilibrium" and see if it can be jus- 
tified under the conditions of (I) and (II) above. We note 
that the collision process involves many time scales, and 
there are at least six of them relevant for our present con- 
sideration: 1. The time scale associated with the infall 
velocity: ti = R/Vi, R=the radius of the NS, Vi— (infall 
velocity at the point of contact). 2. The time scale as- 
sociated with the local sound velocity: t s = R/V s , V s = 
(sound velocity). 3. The time scale associated with the 
velocity of the shock (velocity in the rest frame of fluid) : 
t s h = R/V s h, V s h = (shock velocity). 4. The time scale 
for the merged object to thermalize, in the sense of being 
describable by one single EOS (same K everywhere): t e . 
5. The time scale of neutrino cooling t n . 6. The time 
scale of the gravitational collapse t c . 

Some comments of these time scales are in order. We 
focus on the case of two 1.4 M© NSs. We model them 
with a polytropic EOS P = Kp T with a polytropic index 
of r = 2. The initial K value of the two stars is taken to 
be 1.16 x 10 5 ~~ i 7 . (Maximum stable mass of these values 
of K and T is 1.46 Mq.) We note that the argument in 
is applicable to all polytropic models. 



For this model, Vi is (somewhat larger than) the 
Newtonian value ~ 0.28c, as can be estimated by 
y/GM/(2R); the diameter of the NSs is about 2R = 
26km (the isotropic coordinate radius of this NS is 9.3km, 
the proper radius is R = 13km). Hence the time scale 
associated with the infall velocity ti is about (smaller 
than) 0.16 ms. To estimate the second time scale t s , 
note that the sound velocity V s depends strongly on the 
dynamical process and the region under consideration. 
For the model mentioned above, the initial central rest 
mass density of the NSs is about 1.5 x 10 15 <?/cm 3 ; V s 
there is about 0.5 c. With the density elsewhere initially 
lower than this value, but higher in some period in the 
central region of the collision, V s varies but is roughly 
0.5c. Thus, t s is roughly 0.1ms. To estimate the third 
time scale t s h requires an estimation of the velocity of 
the shock V s h produced in the collision. The locally mea- 
sured proper velocity of the shock V s h is higher than, but 
of the same order of magnitude of, the sound speed V s at 
a fraction of c in the head-on collision case. Hence t s h is 
also of order 0.1ms. These three time scales determine 
the time scale 4 which is central to our discussion. In 
near static situation, or when the bulk velocity of matter 
is small (Vi « V s and Vi « V s h), t e can be taken to 
be a few times t s or t a h- On the other hand, its value 
in a highly dynamic situation with Vi comparable to V s 
and V s h is an important issue to be discussed below. The 
fifth time scale t n governs the final settling down of the 
merged object after t e . t n is of the order of seconds, or- 
ders of magnitude longer than the first four time scales. 
The gravitational collapse time scale t c is in turn con- 
trolled by these time scales 1-5. It can be as short as 
ti, or as long as t n . For the collision of two 1.4 M Q 
NSs, the merged object would have to collapse after t c , 
if not before, for most of the reasonable EOS. We call 
collapse that occurs on the first four time scales prompt 
collapse, and collapse that occurs on a longer time scale, 
like t n , delayed collapse. For more general coalescence 
processes, there can be other time scales involved, e.g., 
the time scale of angular momentum transfer t a , and the 
time scale of gravitational wave emission t g . However, 
for the case of head-on collision with the stars falling in 
from infinity, we expect strong shock heating causing t n 
to be shorter than t g . We do not have to consider t a and 
t g in our present consideration. 

In Shapiro's argument, the time scale 4, t e , is implic- 
itly taken to be the shortest time scale in the problem, 
so that the system can be described by a single EOS at 
any instant in the collision process. The above discussion 
suggests that this may not be true for the two 1.4 M Q 
NS collision case. Indeed, the relations between the time 
scales 1, 2 and 3 strongly affect t e . With ti comparable 
to t s and t s h, dynamic effects are important, and t e can 
be longer than ti. In particular, with matter falling in 
at high speed along the axis of the collision, the speed of 
the shock wave in that direction would be significantly 



reduced, until after ti, delaying "thermalization" of the 
coalescing objects. For situations like this, arguments 
based on a uniform EOS throughout the coalescing ob- 
ject cannot be justified. Indeed, when the infalling time 
scale ti is comparable to the other time scales in the pro- 
cess, it could happen that even if a hydrostatic stable 
equilibrium configuration exist, the dynamics of the sys- 
tem might not lead to that configuration and the time 
scale of collapse could be as short as ti. 

Another way of looking at the problem is to imagine we 
tie the two stars on strings and lower them towards one 
another in a quasi-stationary fashion while depositing the 
potential energy extracted back to the two stars. For this 
case Shapiro's argument would be applicable. However, 
for a NS collision with the time scales discussed above, 
one would have to examine the dynamics of the infall to 
determine whether a prompt or delayed collapse would 
occur. In short, as both a thermally supported merged 
object and a black hole can have the same rest mass and 
total energy, arguments based solely on conservation of 
mass and energy without taking the dynamics into con- 
sideration cannot rule out one outcome from the other. 

We note that the above time scale considerations sug- 
gests that whether it is a delayed or prompt collapse in 
head-on collision can depend on the initial NS's configu- 
ration. It does not imply prompt collapse by itself. To 
demonstrate that a prompt collapse results, one has to 
perform a fully relativistic simulation. 

Our NCSA/Potsdam/Wash U collaboration is devel- 
oping a multi-purpose 3D numerical code, "Cactus" , for 
relativistic astrophysics and gravitational wave astron- 
omy. This code contains the Einstein equations coupled 
to the general relativistic hydrodynamic equations. For 
a description of various aspects of the code and the NS 
grand challenge project based on it see H . Testbeds and 
methods for evolving neutron stars have been given in 
0, and will not be repeated in this paper. While this 
multi-purpose code is still under development for various 
capabilities in treating a broad class of astrophysical sce- 
narios, in this paper we focus on the results obtained by 
applying this code to the head-on collision problem. 

Simulation results. We show the M = 1.4 Mq 
head-on collision case. The stars are modeled as given 
above. We put the two TOV solutions at a proper dis- 
tance of d = 44 km apart (slightly more than 3 R sepa- 
ration) along the z-axis, and boost them towards one an- 
other at the speed (as measured at infinity) of yj GM / d 
(the Newtonian infall velocity). The metric and extrin- 
sic curvature of the two boosted TOV solutions are su- 
perimposed by (i) adding the off-diagonal components of 
the metric, (ii) adding the diagonal components of the 
metric and subtracting 1, and (iii) adding the compo- 
nents of the extrinsic curvature. The resulting matter 
distribution, momentum distributions, conformal part of 
the metric, and transverse traceless part of the extrin- 
sic curvature are then used as input to York's procedure 



H for determining the initial data, in maximal slicing. 
With this setup, the initial data satisfies the complete 
set of Hamiltonian and momentum constraints to high 
accuracy (terms in the constraints cancel to 10 -6 ), and 
physically represent two NSs in head-on collision falling 
in from infinity, at least up to the Newtonian order. (For 
initial data setup on the PIN level, see j6|). 

The initial data is then evolved with the numerical 
methods described in Q]. Various singularity avoiding 
slicings been used and tested against one another (max- 
imal and 1 + log slicings most extensively), yielding ba- 
sically the same results. The simulations have been car- 
ried out with resolutions ranging from Ax = 1.48 km to 
0.246 km (13 to 76 grid points across each NS, with 32 3 
to 192 3 ) for convergence and accuracy analysis. 




Lapse vs. z (km) vs. z (km) 

FIG. 1. a) The lapse (a) and b) g zz along the z-axis are 
displayed at various times. This simulation used 192 3 grid 
points, with Ax — 0.246 km. 



In Fig. la we show the collapse of the lapse along the 
x = y = line from t = ms to t = 0.31 ms at intervals 
of 0.044 ms. (With the reflection symmetry across the 
z = plane and the axisymmetry of the head-on collision, 
we only need to evolve the first octant.) At t = 0.31 ms 
the lapse has collapsed significantly. 

In Fig. lb we show the evolution of g zz along the z 
direction (x = y = 0) . We see the familiar grid stretching 
effect associated with evolving a black hole. 
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FIG. 2. The evolution of the lapse (a), the rest mass 
density (p), and the pressure (P) in the region centered at 
the point x = y = z = 0to the time t = 0.31 ms. 

Fig. 2 shows the time development of the lapse, the 
(proper) rest mass density p and the pressure P at the ori- 



gin, scaled by the critical secular stability values p cr iticai 
and Pcritical, the values beyond which a static TOV so- 
lution is unstable to collapse for the given polytropic co- 
efficient K and index T(= 2). We note that the effective 
K{= P/p 2 ) is time dependent due to shock heating. At 
coordinate time t = 0.26 ms we see that both p and P 
surpass p cn u ca i and P cn ucai, indicating a collapse. 
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FIG. 3. The position of the AH at different resolutions 
and outer boundary locations, all at t = 0.31 ms. 
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FIG. 4. Contour lines of the log of the gradient of the rest 
mass density log ( y^V*(p)Vi(p) J , showing a shock front at 
coordinate radius ~ 4.5 km. 

In Fig. 3 we show the position of the apparent horizon 
(AH). To confirm the location of the AH, convergence 
tests both in terms of resolution and in terms of location 
of the computational boundary have been carried out. 
(For a discussion of the AH finder, see 0]). We have also 
explicitly determined trapped surfaces bounded by the 
AH for the positive confirmation of a collapsed region. 
In Fig. 3, the solid and long dashed lines correspond to 
the AH locations at resolutions of Ax = 0.492 km and 
Aa; = 0.246 km, while the dotted line corresponds to 
Aa; = 0.492 km but with the outer boundary two times 
further out. Although the coordinate position of the AH 
is substantially elongated in the z direction, the AH is 
actually quite spherical. The proper circumference on 
the x-y plane (equatorial) is close to the circumference on 
the x-z plane (polar), with the latter being 52.9±1.9 km. 



For comparison, 4itMah is 52.9 ± 2.1 km, where Mah 
is the mass of the AH (we note that a substantial part 
of the matter in the system is enclosed within the AH) . 
Analysis of this in relation to the hoop conjecture will be 
given elsewhere. 

Fig. 4 shows the contour lines in the y — plane 
of the log of the gradient of the rest mass density 
log ( -\/V*(p)Vi(pH at time t — 0.31 ms. We see a sharp 
peak at a coordinate radius of ~ 4.5 km. The sharp 
change in rest mass density indicates a shock, stronger 
in the infalling direction (z), while weaker near the equa- 
torial plane. The shock is moderately relativistic with 
Lorentz factor of about 1.2. The shock is well captured 
in this 192 3 run with high resolution shock capturing 
(HRSC) GR-hydro treatment. Comparing to Fig. 3, we 
see that the shock front is inside the AH in all direc- 
tion at this time, although it is still moving outward in 
coordinate location. 
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FIG. 5. The evolutions of the L2 norms of (a) the Hamil- 
tonian constraint, and (b) the z-component of the momentum 
constraint. 

In Fig. 5 we show the convergence of the Hamilto- 
nian and the z-momentum constraints for a measure of 
the accuracy of the simulation. The evolution of the L2 
norms (integrated squared) of the constraints are scaled 
by the maximum of the matter terms in the constraints 
(16irp ADM and 8713* respectively). The solid, dotted, 
and dashed lines represent the constraints at resolutions 
Ax = 1.48 km, 0.492 km and 0.246 km respectively. 
These long time scale convergence tests indicate that our 
numerical evolution is stable and convergent for the time 
scale of our present problem. Towards the end we see 
that the error is increasing rapidly; an examination of 
the spatial distribution of the constraint violations shows 
that the error is due to the familiar problem of resolv- 
ing the "grid stretching" peaks of the black hole metric 
(cf. Fig. 2). An extensive convergence analysis of many 
of the variables involved in the simulation has been car- 
ried out and will be presented in a follow up paper. We 
have also performed simulations with the initial boost ve- 
locity increased by 10% (generating more shock heating) 
and confirmed that our results are not sensitive to the 
initial velocity. 

With these results we conclude that prompt collapse 
of the merged object formed in head-on collision infalling 



from infinity is possible, under the same conditions as in 
Shapiro's conjecture. 

We have also carried out simulations of head-on colli- 
sions of lower mass NSs and have seen cases in which the 
shocks propagate to cover the whole star and no AH is 
found, indicating that the collapse would be delayed un- 
til radiative cooling. A detailed analysis of the transition 
point between prompt and delay collapse is computation- 
ally expensive with our 3D code used to carry out the 
present analysis. A 2D version of the present treatment 
is being developed with this specific application in mind. 

Conclusions. We pointed out that there is an as- 
sumption in Shapiro's conjecture, namely, the head-on 
collision process is in quasi-equilibrium (in the sense of 
(A) and (B) above). We showed that this may not be 
true for the collision of two 1.4 Mq NS's. We substanti- 
ated our argument with a simulation solving the full set 
of the coupled Einstein and general relativistic hydrody- 
namic equations. We confirmed the prompt formation of 
a black hole in the infalling time scale ij with an apparent 
horizon found 0.16 ms after the point of contact. 

In this paper we concentrate on the head-on collision 
process under the same conditions as in Shapiro's conjec- 
ture. As the time scale argument given above is rather 
general, and in particular does not depend on the poly- 
tropc EOS, we expect the same argument to be applica- 
ble to more general situations. An investigation of the 
prompt vs. dlayed collapse problem of head-on collisions 
with realistic EOSs, more realistic initial conditions (ini- 
tial data setup with Post-Newtonian formulation), and 
with a determination of the critical point between delayed 
vs. prompt collapse will be given in follow up papers. 
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